Viscoelastic model for the dynamic structure factors of binary systems. 
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This paper presents the viscoelastic model for the Ashcroft-Langreth dynamic structure factors 
of liquid binary mixtures. We also provide expressions for the Bhatia- Thornton dynamic structure 
factors and, within these expressions, show how the model reproduces both the dynamic and the 
self-dynamic structure factors corresponding to a one-component system in the appropriate limits 
(pseudobinary system or zero concentration of one component). In particular we analyze the be- 
havior of the concentration-concentration dynamic structure factor and longitudinal current, and 
their corresponding counterparts in the one-component limit, namely, the self dynamic structure 
factor and self longitudinal current. The results for several lithium alloys with different ordering 
tendencies are compared with computer simulations data, leading to a good qualitative agreement, 
and showing the natural appearance in the model of the fast sound phenomenon. 



I. INTRODUCTION. 

The development of inelastic neutron scattering (INS) 
techniques opened up, around 40 years ago, the exper- 
imental study of the dynamic properties of several con- 
densed matter systems, in particular of liquids. In princi- 
ple the total scattered intensity in an INS experiment in- 
cludes both incoherent and coherent contributions, which 
are related respectively to the self-dynamic structure fac- 
tors and the dynamic structure factors. A clearcut sep- 
aration of both contributions is not always possible and 
in the analysis of the raw data it is useful to have simple 
models for the dynamic and/or self-dynamic structure 
factors in order to achieve such a separation through a 
numerical fitting procedure, and perform a proper inter- 
pretation of the experimental data. Even in those cases 
where there is coherent scattering only, it may happen 
that the particular behavior of the dynamic structure 
factor as a function of frequency obscures the analysis, 
for instance when no clear side peaks appear; in this case 
again the availability of models for fitting helps in the in- 
terpretation of the mechanisms controlling the behavior 
of the dynamic properties of the system. Similar prob- 
lems are encountered when the dynamic properties of liq- 
uid systems are studied by either inelastic X-ray scatter- 
ing (IXS) or molecular dynamics (MD) simulations. Even 
though IXS produces coherent scattering only, and MD 
provides a very detailed information of the properties of 
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interest, nevertheless the interpretation of the numerical 
data obtained is greatly aided if theoretical models are 
available. 

In this respect, and for pure systems (one-component 
systems), a prominent role has been played by the so 
called viscoelastic model, introduced by Lovesey, which 
basically describes the dynamic structure factor as a 
sum of three Lorentzian functions of frequency, one of 
them representing particle diffusion and the other two 
describing damped propagation of collective excitations. 
This model applies for intermediate wave vectors, k, in- 
between those corresponding to a hydrodynamic behav- 
ior (low fc, where the hydrodynamic model is applicable) 
and those corresponding to an ideal gas behavior (large 
k, where the free-particle model is correct). A similar 
expression is also available for the self-dynamic structure 
factor, but its use has been much more scarce in the liter- 
ature, although as we shall show below, the viscoelastic 
model for the self-dynamic structure factor in fact has 
very interesting properties that other models lack. 

The case of mixtures is more complicated. Although 
both the hydrodynamic and the free-particle models are 
readily extended to mixtures, there is no well-behaved 
model so far to describe the intermediate k range. A pre- 
vious attempt to extend the viscoelastic model to liquid 
mixtures , failed because some errors in the derivation 
made it incorrect, and therefore inapplicable. In partic- 
ular the model did not recover the one-component case 
for pseudobinary systems, i.e. those systems which are 
in fact one-component, but with some particles labeled 
different from others (some of them are named type 1 
and the rest type 2). 

In this paper we extend the viscoelastic model to mix- 
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tures, and in particular to binary mixtures, in a con- 
sistent way which reduces to the correct one-component 
hmit in the appropriate cases (pseudobinary mixture or 
zero concentration of one component). Moreover we give 
expressions for the so called Bhatia-Thornton dynamic 
structure factors, which are very useful when discussing 
ordering properties of binary systems. From these ex- 
pressions it is easy to show the reduction to the one- 
component case, leading not only to the dynamic struc- 
ture factor, but also to the self-dynamic structure factor 
of pure systems. 

We compare the results of the model with MD data for 
three different type of systems: the first is a pseudobi- 
nary alloy, namely, pure liquid Li, the second corresponds 
to Li-Mg, which is a typical quasi ideal system, and fi- 
nally the third one corresponds to Li4Pb, which is an 
archetypical case for a class of compound forming alloys. 
The case of systems with tendency to phase separation 
has already been considered before in the study of 
liquid Lio.6iNao.39, which is again a typical phase sepa- 
rating mixture. Although apparently trivial, the study of 
the pseudobinary case leads to interesting conclusions re- 
garding the behavior of the self-dynamic structure factor, 
which vindicate the use of the viscoelastic model for the 
self-dynamic structure factor of one-component systems. 



II. FORMALISM. 

A. One-component system. 

Here we merely recall the expressions for the dynamic 
properties we are interested in for one-component sys- 
tems. The basic magnitude to be considered is the in- 
termediate scattering function (ISF), F{k,t), which de- 
scribes the collective dynamic behavior of the system and 
is defined as F{k,t) = (pg(i)p_j(0)), where Pj^{t) is the 
microscopic number density. 
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initial values Fs{k, 0) = 1, Fd{k, 0) = S{k) - 1. The SISF 
is of interest by itself, since it is the basic magnitude 
that describes the one-particle dynamic behavior of the 
system. 

The dynamic structure factor (DSF), S{k,uj), and the 
self dynamic structure factor (SDSF), Ss{k,uj), are ob- 
tained from the corresponding intermediate scattering 
functions by passing to the Fourier domain: 
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S{k, uj) = — dte"^* F{k, t)^- Re F{k, z = 

(3) 

where Re stands for the real part and F{k, z) is the 
Laplace transform of F{k, t), i.e. 



F{k,z) 



dt e"^* F(k,t) 



(4) 



The memory functions of the ISF, Af(fc, i), and of the 
SIFS, Ms{k,t), are then introduced through the general- 
ized Langevin equations, which read in the time domain 
and in the Laplace domain: 



dr M{k,T) F{k,t-T) 



F{k,z)= z + M{k,z) F{k,t = 0) 



(5) 
(6) 



with equivalent equations for the self counterparts. The 
higher order memory functions are introduced exactly 
in the same way: the second order memory functions 
N{k, t) and A^s(^, t) are the memory functions of M{k, t) 
and Ms{k,t) respectively, while the third order memory 
functions, K{k, t) and Ks{k, t) are the memory functions 
of N{k,t) and Ns{k,t) respectively. The initial values of 
the memory functions which appear in the Laplace for- 
mulation of the generalized Langevin equation are easily 
determined in terms of the nth frequency moments of the 
DSF and SDFS, w"(fc) = /.!^ S{k,uj) duj, 



of the system composed of N particles at positions re{t), 
which are enclosed in a volume V, so that the ionic num- 
ber density is p = N/V . The initial value of the ISF is 
the static structure factor, S{k), which is directly related 
to the pair distribution function, g{r): 

F{k,0) = S{k) ^ 1+p J dr{g{r) - l)exp[-i k-r\ . (1) 

A similar separation into two terms is also possible for 
all times, 



F{k,t)^ F,{k,t)+Fd{k,t) 



(2) 



where we have introduced the self intermediate scattering 
function (SISF), Fs{k,t), and the distinct intermediate 
scattering fimction, Fd{k,t), which obviously have the 



M{k,t^O) = uj^{k) uj"(k)-^ (7) 
N{k,t = 0) = Lu%k)u;^{k)-^ - iu^{k) Lu^ik)-^ (8) 

with again equivalent expressions for the self counter- 
parts. Moreover these frequency moments up to the 
fourth can be determined from the knowledge of only the 
temperature T, the atomic mass, the interatomic pair 
potential, (pir), and the static structure, i.e., g{r) and 
Sik). 

A useful feature of a memory function is that it de- 
cays in time more rapidly than the function from which 
it originates. Based on this, it seems plausible that at 
some level of the hierarchy of memory functions an ap- 
proximation where the memory function is just a Dirac 
delta function at t = should be a good ansatz. In the 
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viscoelastic model the level at which this ansatz is taken 
is the third, i.e., it is assumed that 



K{k,t) = K{k)6{t) 



K{k,z) = K{k) 
Ks{k,z)^ K,{k) (9) 



Explicit expressions of K{k) and Ks{k) in terms of 
the same magnitudes as the frequency moments can be 
obtained 0, 3 by imposing that in the free particle limit 
(fc — > cxd) the theory recovers the known exact S{k,uj — 
0) = Ss{k,uj — 0). Further details are given in the next 
section. 

Introducing in eq. ^ the higher order memory func- 
tions up to the third it is found that 



F(k,z) 



F{k,t^O) _P2(k,z) 
M{k,t^O) ~ P3{k,z) 



(10) 



N{k,t = 0) 



z + 



z + K{k) 



where Pn{k,z) denotes a real polynomial in z of degree 
n. Making a partial fraction decomposition of this ex- 
presion, and denoting by Zi the roots of Ps{k, z), which, 
by the way, are the same as the roots of z + M{k, z), we 
have 



F{k,z) = 



3 



F{k, t)^^A, exp[zii] 

i=l 



(11) 

where the Zi and the Ai appear either as real or in com- 
plex conjugate pairs. Therefore either the three roots are 
real or one is real and the other two are a complex con- 
jugate pair. In all practical situations the latter is the 
case, and then the roots are interpreted as describing a 
diffusive mode and a pair of damped propagating modes, 
much the same as in the hydrodynamic model, although 
for k values outside this regime. The DSF is then a sum of 
three Lorentzian functions, one centered &t uj — Q which 
corresponds to the real root, and the other two centered 
at the imaginary parts of the complex conjugate pair. 

At small fc, the functional form of the dynamic mag- 
nitudes within the viscoelastic model coincides with that 
of the hydrodynamic model, which is known to be accu- 
rate in this region. Futhermore, the viscoelastic modes 
behave (as a function of k) much in same way as the 
hydrodynamic modes. However, in spite of the previous 
similarities, there are basic differences between the vis- 
coelastic and the hydrodynamic models. They are bet- 
ter understood if the derivation of both models is made 
by an alternative route using the generalized modes ap- 
proach 0. Here one considers the equations of motion 
of several dynamic magnitudes, namely temperature (or 
energy) fluctuations, density fluctuations, current fluc- 
tuations and stress tensor fluctuations. The viscoelastic 
model follows from considering couplings among the last 
three variables and ignoring their coupling with tempera- 
ture fluctuations. The hydrodynamic model follows from 



considering couplings between the first three magnitudes 
and ignoring their coupling to the stress tensor fluctu- 
ations. A parameter quantifying the coupling between 
density and temperature fluctuations is the specific heat 
ratio 7 — Cp/Cv If 7 ^ 1 the coupling is weak and the 
viscoelastic model is expected to be accurate; otherwise 
the model is expected to fail especially at small k. 

Within the viscoelastic model, the analytical structure 
of the self functions, i.e. Fs{k, z), the SISF and the SDSF, 
is the same as that of the equivalent collective counter- 
parts. However, the amplitude of the real coefhcient as- 
sociated to the real root is usually much greater than the 
amplitudes of the complex terms and therefore the asso- 
ciated propagating modes in the SDSF have often been 
neglected when the viscoelastic model has been applied 
to study the dynamic properties of one-component liq- 
uids. It may be argued whether these modes are real or a 
mere artifact of the model because the SDSF describe the 
one-particle behavior of the system. However, the close 
connection between collective and single particle prop- 
erties in dense systems could induce the appearance of 
these modes. In any case, we stress that, to our knowl- 
edge, no detailed study of this possibihty has yet been 
carried out. 



B. Binary system: Ashcroft-Langreth partials. 

The generalization of the foregoing formalism to binary 
systems is straightforward. We consider N particles in a 
volume V , therefore with number density p = N /V , com- 
posed of two species (i = 1, 2) which are characterized by 
the number of particles Ni, or concentration Xi = Ni/N, 
the atomic masses m^, and ionic partial number densi- 
ties Pi — Xip. The interaction between particles of type 
i and j is described by effective interatomic pair poten- 
tials <pij{r) and the static structure by the partial pair 
distribution hmctions gij (r) . 

We start from the Fourier transform of the partial mi- 
croscopic number densities, defined as 



E 

3 1=\ 



exp i k-r-^^it) {j = 1,2) 

(12) 

from which the partial intermediate scattering functions, 
Fij{k,t) are obtained as 



F,ik,t)^<f^(t)p^''(0)>= F(fc,t) 



0), 



(13) 



where the last equation defines the 2x2 matrix F(fc, t). 
The initial value of this matrix reduces to the matrix of 
Ashcroft-Langreth (AL) partial structure factors. 



F(fc,t = 0) = S(fc) 



(14) 
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where 

[S(^)],, = S,,ik) 

= Sij + {xiXjY^'^p J df {gij{r) - l)exp[-z k ■ f] 

(15) 

As in the case of one-component systems, a similar sep- 
aration into two terms is possible for all times, defining 
the self and distinct parts of the partial ISFs: 

F(fc,t) = Fs(fc,i) +Fd(fc,t) 
Fy(fc,t) = S^jF,^j{k,t) + {x,XjY^^Fd,^,{k,t) (16) 

Equations lO-® remain formally the same, although 
now all the magnitudes in the equations are matrices, 
the products are to be understood as matrix products 
and the exponent -1 means matrix inverse. For instance, 
the memory function matrix obeys in the time and the 
Laplace domain the equations 



—F{k,t)^- / dTM{k,T)F{k,t- t) (17) 



F{k,z)^ zI + M{k,z) F(fc,< = 0) (18) 

where I is the unit matrix. 

The matrices of frequency moments are given by |^ 



w°(fc) -5,y(fc) 



\k) 



^ksT 



(19) 



UJ^{k) 



Sij k 



knT 



mifTij 

-{xiXj^ ^ k 

,2 ksT 



ik'^ksT + '^xep / dr gui; 



WLiraj 



p df Qij (r) cos(fcf^ 



rmmj 



ik^ksT + ^ XiXiiik) 



(x,a:,y/2xL;y(fc)} 



(20) 



where T is the temperature, ks is Boltzmann's constant, 
z denotes the direction of k and the last equality defines 
the functions Xij and XL-ij- 

The hierarchy of memory function matrices is again 
truncated at the third level, setting 



K{k,t) = K{k)S{t) 



K{k,z)=K{k) (21) 



so that explicitly we have for the ISF matrix and the first 
and second order memory function matrices the relation 
(|TH|l and 



M{k,z)= zI + N(fc,z) M{k,t^O) 



(22) 



N{k,z)= zI + K{k) N(fc,i = 0) 



(23) 



with the initial values of the memory function matrices 
given in terms of the frequency moments matrices as in 
eqs. Q and ©. In order to provide explicit expressions 
for the matrix elements of K(fc) we follow Lovesey's ar- 
guments 13,0. Setting z = in eqs. and we 
have 

M(fc,z = 0) = N(fc,z = 0)"^ M(fc,i = 0) 

= K(fc)-^ N(fc,i = 0)1 ^M(fc,t = 0) 
^ N{k,t = 0)-^ K{k)M{k,t^O) (24) 

and therefore 

K(fc) = N(fc, t = 0) M(fc, z = 0) M(fc, t = 0)-^ (25) 

The value of M(fc, z = 0) = dt M(fc, t) is then as- 
sumed to be rather insensitive to the detailed shape of 
the Mij{k, t) and therefore, a reasonable estimate can be 
obtained from an approximation that fulfills its correct 
short time behavior obtained by a simple Taylor expan- 
sion, namely. 



M(yfc,t) = 



I N(k,t 

2 ^ ' 



0) 



M{k,t = 0) 



(26) 

The specific approximation we make for M(fc,t) to per- 
form the integral is 



M(fc, t) « exp 



-|-N(fc,< = 0) 



M{k,t = 0) . (27) 



It is precisely at this point where the previous attempt to 
generalize the viscoelastic model [J failed, since a sim- 
ilar approximation was made at the level of the matrix 
elements Mij{k,t) only, and not for the matrix M(fc, t) 
itself, which is obviously wrong because we are dealing 
with matrix products. 

Within this approximation M(fc, z = 0) would be given 
by the time integral of the exponential function times the 
matrix M(fc,t — 0). Trying to correct the inaccuracies 
that might have been introduced by this approximation, 
this value is premultiplied by a matrix H of constants to 
be determined later. In this way an explicit expression 
for K(fc) is found in terms of the initial values of the 
second order memory functions: 

- - r ^2 ■ 

K(fc) = N(fc,i = 0) H / exp N(fc,t = 0) 

"'0 [2 

(28) 

The determination of H is carried out by imposing that 
in the free particle limit (fc — > oo) the dynamic structure 
factor matrix evaluated at zero frequency recovers the 
exact free gas result, i.e.. 



[S(fc,c^ = 0)]y =(5,. 



/ 1 



1/2 



V27rfc2 ksT J 



(29) 
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leading to the result H = [2^/2 /tt)!. In the Appendix we 
give details of this derivation and explicit expressions for 
the matrix elements of K(fc). 

Turning now to the analytic behavior of F(fc,z) we 
see in eq. (I18|l that it is determined by the inverse of the 
matrix [zI+M(fc, z)], namely the transpose of the adjoint 
divided by the determinant. Therefore, all the Fij{k,z) 
have the same analytic behavior, which is determined by 
the solutions of the equation 



det 



zI + M{k,z] 







(30) 



As in the one-component case, this determinant is a 
real rational function, whose form is obtained by writing 
down explicitly the equations for the memory function 
matrices up to the third, leading now to an expression 
of the type Pe{k, z)/Pi{k, z). Consequently there will be 
six roots, which in principle may appear (together with 
the corresponding amplitude matrices) as 6 real values, 
4 real values and a complex conjugate pair, 2 real values 
and 2 complex conjugate pairs, or 3 complex conjugate 
pairs. Therefore we have 



F(fc,z) 



6 

i=l 



1 



- Ai F(fc, t) = ^ Ai exp[zii] 



1=1 



(31) 

where the Ai are the fc-dependent 2x2 amplitude ma- 
trices. In the actual calculations made for liquid Li-Na 
0, as well as those carried out in this paper, we found 
for all k two real roots and two pairs of complex con- 
jugate roots, at variance with some other calculations 
where one of the complex conjugate pair transforms into 
two real roots below a certain small k value 0, |^ . AH 
roots have negative real parts 0, so we can rewrite the 
previous equation as 



F^,{k,z) = 



A* 



B* 



- z(i)* 



- 7P) 



- z(2)* 



0-hz(3) z + z(4)) 



(32) 



where the fc-dependent coefficients Aij and Bij are com- 
plex, while Cij and Dij are real. The complex roots, 
z(j)(fc) = -r(j)(fc) ± iJ^\k), {j = 1,2) describe propa- 
gation, where the real part r(^-'(fc), represents damping 

whereas the imaginary part Lui''\k), represents the prop- 
agation. A set of two complex conjugate roots repre- 
sents propagation in opposite directions. The real roots, 
—z^^\k), {j = 3,4) stand for purely diffusive processes. 
Out of the six roots, only four go to zero when fc ^ 0; 
this which coincides with the behavior of the four hydro- 
dynamic roots predicted by the hydrodynamic equations 
for binary mixtures. Again we stress that although the 
small fc behavior of the four viscoelastic modes is similar 
to the hydrodynamic ones, there are quantitative differ- 
ences between them based on the same reason as in the 
one-component case, namely the neglect of coupling with 



the energy fluctuations. As in one-component liquids, the 
viscoelastic model is expected to work for systems where 
this coupling is weak, i.e., those with specific heat ratio 
7 ~ 1. The other two viscoelastic roots have a finite 
value when k —^ 0, while the corresponding amplitudes 
vanish in this limit; this is the typical behavior of kinetic 
modes. The six roots lead to a Sij{k,uj) given as a sum 
of six Lorentzians, two centered at w = and the other 
four at w = ±uji^\k), (j = 1, 2). 



C. Binary system: Bhatia-Thornton partials. 

An alternative description of the structure of binary 
alloys is provided by the so called Bhatia-Thornton 
(BT) functions. The BT static partial structure factors, 
namely, number- number SNNik), number-concentration 
SNcik), and concentration-concentration Sccik) partial 
structure factors describe the correlations among fluctu- 
ations in number density (irrespective of chemical compo- 
sition) and concentration density, and are linear combi- 
nations of the AL partial structure factors Sij (fc) . These 
relations are most compactly written if one defines the 
matrix S{k) of modified BT partial structure factors in 
terms of the matrix S(fc) of AL partial structure factors: 



where 

m 

and 



5 = X S X 



SNNik) SNcik) /{XIX2)^/^ 

ScNik)/ixiX2)'^/'^ Sccik) /ixiX2) 



X = x^ = 



(33) 



(34) 



Xl 

x^ 



X2 
X'l 



All other BT magnitudes (intermediate scattering 
functions, first, second and third memory functions, fre- 
quency moments, etc.) are defined in the same way, pre- 
and post- multiplying the corresponding matrix of AL 
magnitudes by the matrix X; for instance, the BT par- 
tial ISPs are defined by the matrix J^(fc,t), 

Tik, t) = X F(fc, i) X = Tsik, t) + Tdik, t) (35) 
which leads to (dropping the arguments of the functions) 

FnN = {xiFs^l + X2Fs^2} 

+ {x\Fd,ii + 2xiX2Fd^i2 + xlFd^22] (36) 



F, 



Nc 



{ixiX2)'/'iFs,l~F.s,2)} 



ix\X2)^l'^ 

+ {(xixa)^/' \xiiFisi - Fis2) + X2iFi^^2 - i^d,22)]} 

(37) 
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and 



^ cc 

X1X2 



+ {xiX2[Fd,ii-2Fd,i2 + Fd.22)]} (38) 



where the first braces in each equation enclose the BT 
SISFs and the second braces enclose the BT distinct ISFs. 

The viscoelastic model for the BT functions is then 
defined by the relations 

:F{k,z)^\zl + Mik,z)\ V(fc,i = 0)(39) 

-1 



M{k,z)= zI + Af{k,z) 
Af{k,z) = \zI + IC(k) 
/C(fc) = 7V(fc, < = 0) H dt exp 



M{k,t^ 0) (40) 
V(fc,t = 0) (41) 



-Af(k,t^ 0) 



(42) 



T{k,t = Q) = ^^{k)^S{k) (43) 
M{k,t = 0) = n^{k)n°{k)-^ (44) 

A/'(fc, t = o) = n-^{k)n^{k)-^ - n'^{k)n°{k)-^ (45) 

rj"(fc) = X w"(fc) X (46) 



explicitly 



ni 



k^keT 



^Nc = k'^kBT{xiX2Y 



and 



7712 
1 

mi 



1/2 ^ L 



X2_ 
7772 
X2_ 
?77l 
1 
7772 



(47) 



= 2 3fc /i;BT + a;iX2 X12 



K ksl \7nf 'Tio / V'^2 Til 



2 2 \ 

^1 r 1 c\ X\X2 r 1 2^2 r i \ 

— [Xll - XL;ll\ + 2 [Xl2 ~ XL;12\ + — [X22 - XL;22 

771J 77717)7,2 "^2 / 



o4 / \ / \ 2 

"cc / 2;i , a;2 



7)7,S 



j 3PkBT + 








\m-2 


mi J 



Xl2 



+ X1X2 ( [Xll - XL;ll] [Xl2 - XL;12] + ^ [X22 " XL;22 

^Nc _ ( ^ ^ \ il2i^ rr , ^ ( Xll Xl2^ , _ /^Xl2 X22 



k^kBT{xiX2Y/^ 



7)7,1 



3fc ksT + xi — 5- ^ + a;2 



7)7 1 7772 




/ Xll 




\ 777i 


777^7 


/ XL;12 


Xi:22 


\ 777 l 7772 


7772 



xi 2 ^2 2- (48) 

' 7)7f 777 l7)72 ' I ™ ™ ™^ / 



D. Reduction to the one-component case. 

In the case of a pseudobinary mixture, when all the 
particles are of the same type but are labelled different, 
all the partial pair distribution functions reduce to that 
of the real one-component liquid, i.e., 311(7") = g22{i') = 
9i2{r) = g{r). However, the AL partial structure factors, 
do not coincide with that of the real one-component sys- 
tem, S{k), but 

S^j(k)^5,, + {x,x,Y'''[S[k)~l] . 

Explicitly we have 

S'ii(fc) = X2+ xiS{k) 

S22ik) ^ Xi+ X2S{k) 

Si2{k) = -{xiX2y/^ + {xiX2f'S{k) . (49) 



The BT partial structure factors, on the other hand, be- 
have more simply, because we now have SNN{k) = S{k), 
SNc{k) = and 5'cc(fc)/(a;iX2) = 1. 

In many one-component systems (in particular liquid 
metals near the triple point) the value of S{k) for small 
k is rather small, of the order 2 — 3 x 10"^, reflecting 
a small value of the isothermal compressibility. For a 
pseudobinary system, we therefore obtain that, except 
for very dilute mixtures, the small k values of the partial 
static structure factors are dominated by the first terms 
of eqns. (|49|l . which take on significantly higher values 
than that of the structure factor of the one-component 
liquid. 

The situation concerning the ISFs is similar. Accord- 
ing to their definition, see eqn. we have Fs.i(fc, t) = 
Fsa{k,t) = Fs{k,t), which is the SISF of the real 
one-component liquid, and Fd.ii{k,t) — -Frf. 22(^1 = 
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Fd,i2{k,t) = Fci{k,t), which is the distinct ISF of the real 
one-component hquid. Therefore the AL partial ISFs are 
given by 

and explicitly, 

Fn{k,t) = X2Fs{k,t) + xiF{k,t) 
F22{k,t) = xiF,{k,t)+X2F{k,t) 
Fi2ik,t) = -{xiX2y^^F,{k,t) + ix,X2y^^Fik,t) , 

(50) 

which is an ackward relation that induces a behavior of 
the partial ISFs very different from that of the real one- 
component liquid. On the other hand the BT partial 
ISFs become 

FNN{k,t) = {xi+X2)Fs{k,t) + {xi+X2fFd{k,t) 
= F,{k,t) + Fd{k,t)^F{k,t) (51) 
FNc{k,t) = (52) 
F,,{k,t) /{xiX2) = ixi+X2)Fsik,t)^ Fs{k,t) , (53) 

that is, the cross function vanishes, the number-number 
function reduces to the ISF of the real one-component liq- 
uid, and the concentration-concentration function, prop- 
erly normalized, reduces to the SISF of the real one- 
component system. 

Coming to the viscoelastic model, in particular in the 
BT formulation, we note that all the magnitudes are de- 
termined by the frequency moments matrices. In a pseu- 
dobinary system we have mi = m2, all the Xij reduce 
to X which is obtained by equalling gij{r) — g{r) and 
4>ij(r) = 0(r); also by the same procedure, all the XL-ij 
reduce to xl- Then all the frequency moments matrices 
in the BT formulation become diagonal, and therefore 
the same occurs to all the other matrices, so that the 
matrix operations in fact become the same operations on 
the diagonal elements of the matrices, i.e., the NN and 
the cc functions. 

In particular we have 



NN 



^J^[3k^kBT + x(k)-XL{k)] (54) 



4 ^'L^[3k'kBT + x{k)] (55) 



which are respectively the 4th frequency moments of the 
DSF and the SDSF of the one-component system. Like- 
wise 



2 _ fc^fcgT _ 2 



(56) 



which again are the second frequency moments of the 
DSF and the SDSF of the one-component system, and 
as we stated before = S^^ik) — S{k) and il^c — 

Scc{k)/{xiX2) = 1, which are the corresponding zeroth 



frequency moments of the DSF and the SDSF of the one- 
component system, respectively. 

Therefore, the present formulation of the viscoelastic 
model for the collective dynamic properties of binary 
mixtures recovers, in the case of a pseudobinary system, 
not only the viscoelastic model for the collective prop- 
erties of the underlying one-component system, through 
the number-number BT functions, but also the viscoelas- 
tic model for its single-particle properties, through the 
concentration-concentration functions. Moreover, the de- 
terminant that leads to the different modes factorizes into 
two terms, one that includes the modes of the DSF with 
the other term accounting for the modes of the SDSF of 
the one-component system; this implies that both sets of 
modes are decoupled. 



III. RESULTS. 

In this section we study three systems and compare 
the results of the viscoelastic model with those obtained 
by molecular dynamics (MD) simulations. 

The first is a pseudobinary system, in which all parti- 
cles are in fact the same type, but half of them are labeled 
as 1 and the other half as 2 {xi = X2 — 0.5). Specifically 
the system is representative of ^Li at T = 725 K and 
p = 0.042 A~'^, and the effective pair potentials used (all 
the three identical) were obtained using the empty core 
pseudopotential 0, with core radius 1.44 A. Some MD 
results for this system have been reported elsewhere ^3 
and are extended here. They have been obtained using 
668 particles in a cubic box with periodic boundary con- 
ditions. 

The other two cases correspond to two Li-based al- 
loys, Lio.yMgo.a, which is a typical quasi-ideal mixture, 
and Li4Pb, which is a reference case for a class of com- 
pound forming alloys with pre-ionic ordering. The tem- 
peratures and densities are T — 887 K and p = 0.040711 
for Lio.yMgo.s and T = 1085 K and p = 0.04558 
A~'^ for Li4Pb. The interatomic pair potentials were ob- 
tained in the case of Lio.vMgo.s within the neutral pseu- 
doatom method 11] while in the case of Li4Pb were taken 
from Copestake et al j^^l- The collective propeties of 
the Lin.yMgo.a alloy have been studied by Anento et al 
[Tsl Il3 | , and those of Li4Pb by several authors 0, 0| , 
the latter being the first system where the fast sound phe- 
nomenon was observed, which consists in the appearance 
of a peak in the Li-Li dynamic structure factor, absent in 
the Pb-Pb one, with a very high frequency similar to that 
of pure Li. In this paper we include MD results obtained 
for both systems, in the case of Lio.rMgo.s with 570 par- 
ticles and in the case of Li4Pb with 600 particles for most 
k values and with 11000 particles for the smallest ones 
(fc„,in = 0.10 A-1). 

The difference between the AL ISFs of the two type of 
systems is striking: while Li4Pb has Fij(k, t), and in par- 
ticular the Li-Li partial, that decay more or less rapidly 
with time, similar to the behavior of one-component 
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systems, Lio.rMgo.s on the other hand shows Fij{k,t) 
which have much larger magnitudes and decay much 
more slowly, in practice making very difficult the Fourier 
Transformation that leads to the DSFs. 

A proper description of such different behaviors, along 
with the case already studied of phase-separating systems 
like Lio.69Nao.31 5) I would therefore provide reliability to 
the viscoelastic model. 

The static structure and the frequency moments, which 
are a required input of the viscoelastic model for the 
calculation of the dynamic properties, have been com- 
puted using the variational modified hypernetted chain 
(VMHNC) approximation which reproduces rather 
well the simulation results, although small differences can 
and do appear in the structural functions, which show up, 
for instance, in a small mismatch of the initial theoret- 
ical and simulation values of the ISPs. While it would 
be possible to use the simulated Sij{k) instead of the 
VMHNC ones, the calculation by simulation of the fourth 
frequency moments is rather unreliable. 

We recall that the viscoelastic model is expected to be 
accurate for systems with 7 '--^ 1. Generalized hydrody- 
namics calculations of 7 have been performed for liquid 
Li [11 , liquid Pb [H, and liquid Na-K and K-Cs alloys 
|20|. leading in all cases to 7 values not far from 1. These 
results give some support to the application of the vis- 
coelastic model to liquid metals and alloys, in contrast 
to other systems like Lennard Jones liquids 0, where 
7 is larger. 



A. Model 1. Pseudo-binary system. 

Due to the particular concentration xi = X2 — 0.5, an 
additional symmetry appears in the system, that implies 
that 5ii(fc) — 522 (fc), Fii{k,t) = F22{k,t) and so on. 

Figure shows the partial structure factors obtained 
form the simulation and the VMHNC theory, which are 
practically coincident. The first thing to note is that, as 
stated above, the small k values attained by the Sij{k) 
are in magnitude close to 0.5 (the concentrations) which 
is one order of magnitude larger than the structure factor 
of the corresponding one-component system in the same 
region. This will imply a much larger initial value of 
the partial ISFs of the system as compared to the ISF 
of the one-component system. The Fij(k, t) are shown in 
fig. 121 where we have plotted both the MD results and the 
viscoelastic ones. The results for small k are markedly 
different from the F(k,t) of the one-component system, 
and not only in the initial value. The latter shows clear 
oscillations around zero, whereas the partial ISFs of the 
pseudobinary system show a large diffusive component 
which decays very slowly. 

The reason for this behavior can be found in the roots 
and corresponding amplitudes that are obtained within 
the viscoelastic model. The roots are plotted in fig. O 
where we see that one of the propagating roots vanishes 
in the fc = limit, while the other remains non-zero 



in this limit. We recall that in a pseudobinary system 
the determinant that leads to the roots factors into two 
third degree polynomials, which correspond to the modes 
of the DSF and the SDFS of the one-component system 
respectively. The first propagating root corresponds to 
the propagating mode in the DSF of the one-component 
system, while the non-vanishing root is an unexpected 
propagating mode that appears in the SDSF of the one- 
component system. Of the remaining two real roots, the 
largest one in magnitude for small k is the diffusive mode 
of the DSF of the one-component system, whereas the 
smallest one in magnitude corresponds to the diffusive 
mode of the SDSF of the one-component system. It is 
precisely this last root that is responsible for the behav- 
ior of the partial Fij{k,t), since its contribution is the 
most slowly decaying one. Moreover its amplitude is the 
largest for small k, as can be observed in fig.0| In fact we 
have only plotted the amplitudes corresponding to the 11 
partial, since due to the symmetry imposed by the con- 
centrations the amplitudes corresponding to the 22 par- 
tial are the same as the 11 ones, while we also have that 
A12 = Aii,Bi2 = -Bii,Ci2 = Cii,Di2 = -Dii. Note 
that Bii goes to as fc ^ 0, which is another character- 
istic feature of the kinetic modes. It is also interesting 
to observe the similarity between Cn and the structure 
factors. We can therefore conclude that the extremely 
slow decay of Fijlk,t) is a direct consequence of the in- 
fluence of the single-particle dynamics on the behavior 
of the ISFs, and this is due to the very definition of the 
AL partial ISFs. When one comes to the BT partials 
no such problems arise, since there is a complete decou- 
pling of the single particle dynamics, which goes exclu- 
sively into the concentration-concentration partial, and 
the collective dynamics, which is exclusively responsible 
for the number-number partial. Obviously this complete 
decoupling is due to the ideal character of this system. 

The Fourier transforms of the longitudinal current 
correlation functions are directly related to the DSFs, 
through the relation Cijlk,uj) = ^Sij{k,ijj) (and the 
same relation for the BT currents), while they can be 
calculated independently in the MD runs. They are es- 
pecially useful in the cases where a slowly decaying term 
appears in the ISFs that complicates the computation of 
the DSFs, while the multiplication by depletes the low 
frequency modes and enhances the high frequency ones. 
In fig.|31we show the functions Cy(fc,w) and also the BT 
ones. We always find a clear peak and a minimum in 
Ci2(fc,w), while there is a clear maximum in Cii(fc,w) 
at the position of the peak of Ci2{k,uj), and sometimes 
a shoulder at the position of the minimum. These are of 
course the remnants of the two propagating modes. The 
coincidence of the shoulder of Cii(fc,w) with the mini- 
mum of Ci2(fc,w) suggests that the second propagating 
mode on top of being kinetic is optical in character, with 
the two "species" moving in opposite directions. The 
appearance of an optical mode in a pseudobinary liquid 
system resembles a similar effect that appears for crys- 
tals: when a linear chain with equilibrium distance a is 
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considered as a chain with equihbrium distance 2a with 
a basis formed by two "different" atoms separated by 
a distance a jlJl, then the first Brillouin zone becomes 
half the original one, and the original dispersion relation 
becomes folded into the new one leading to the optical 
mode. The BT currents on the other hand always show a 
single peak, in the NN case (which is the current of the 
one-component system) at the positions of the peaks in 
the 11 and 12 functions, and in the cc case (which is the 
self-current of the one-component system) at the position 
of the minimum of the 12 function. 

The appearance of the peak in the self-currents of the 
one-component system is very clear in the MD simula- 
tions, showing that also in the single-particle properties 
of this system there are indeed propagating contribu- 
tions, of kinetic character. These modes had not been 
reported previously, at least to our knowledge, since the 
focus has usually been put on the initial value and the 
width of the SDSF of the systems, analyzing the behav- 
ior of these properties within the different theoretical ap- 
proaches. The reproduction of these propagating modes 
by the viscoelastic model (for the single-particle proper- 
ties) is therefore an interesting property of this theory, 
and warrants a wider application of the model in the 
analysis of the single-particle properties of, at least, this 
type of dense liquids. 



B. Li-based alloys: ideal and compound forming 
mixtures. 

In the following we will denote Li as component 1. The 
partial static structure factors of both liquid alloys are 
shown in fig. Those corresponding to Lio.yMgo.a are 
similar in character to those of the pseudobinary alloy, 
except for the difference in concentration, so the values 
of Sij (k — 0) are larger than that of a one-component 
system. On the other hand the results for Li4Pb are very 
different, in particular the values of Sij(k — 0) are all the 
three similar in magnitude to the case of one-component 
system; moreover the dip in S'i2(fc) and the coincident 
positions of the maximum of 6*11 (fc) and the minimum of 
S22{k) is indicative of a kind of ionic ordering in the melt 
[2 2.] . The corresponding ISFs are also strikingly different 
for the two systems, especially for small fc, as observed 
in figs.[3and|Hl As in the pseudobinary system, the ISFs 
of Lio.yMgo.a have a very slowly decaying diffusive term 
for small k. In the case of Li4Pb, the slowly decaying 
term is practically absent in Fii(k,t) and while there is 
indeed such a term in Fi2{k, t) and especially in F22{k, t), 
it decays much more rapidly than in Lio.yMgo.a (observe 
the y-axis scales in figs. H and IHJ even though the Pb 
ionic mass is much larger than the Mg one. 

Again the reason for such behaviors can be traced back 
to the roots and the amplitudes of the modes associated 
to the three partials. The roots are shown in fig. We 
see that, due to the mass difference, the frequency ioi^^ 
and damping F*^^) of the low frequency root are larger for 



Lio.rMgo.s than for Li4Pb, while the high frequency root 
and the larger real root have quite similar magnitude and 
variation with k for the two systems, being much more 
influenced by the masses of the components and the ther- 
modynamic state than by the particular structure of the 
alloy. The smallest real root, z'-'^\ also behaves different 
for Lio.yMgo.s and Li4Pb, staying much closer to F*^-'^) for 
the latter. For small fc, in particular, z'^' is significantly 
larger for Li4Pb than for Lio.yMgo.a , which explains the 
faster decay of Fi2{k,t) and F22{k,t) in Li4Pb. 

Figs . 1 1 OlT^ show the amplitudes for Fij{k,t). The first 
noticeable feature is that the F22{k, t) functions are com- 
pletely dominated for all k by the smallest real root, 
and to a lesser extent by the low frequency root, the 
other two amplitudes being smaller, except in the case 
of Lio.yMgo.s in the region of the main peak of the cor- 
responding structure factor and for small fc, where the 
amplitude corresponding to the larger real root is also sig- 
nificant. The case of the Fi2{k, t) functions is rather sim- 
ilar, with the difference that for small k in Lio.yMgo.s the 
three amplitudes corresponding to the two propagating 
roots and the larger real root are of similar magnitude. 
Note again the striking similarity between the amplitudes 
of the smallest real root and the structure factors, with 
D12 and D22 following respectively 6*12 (fc) and 6*22 (fc)- 
The different behavior of these structure factors for the 
two systems implies now that the amplitude of the diffu- 
sive mode is much smaller for Li4Pb than for Lio.yMgo.s. 
As for the Fii(fc,t) functions, the structure of the am- 
plitudes is richer. For medium and large k the dominant 
term is z^^\ the larger real root, with its amplitude Cn 
following 5ii(fc). For small k, on the other hand, the 
situation is different; Dn is dominant for Lio.rMgo.s , 
while for Li4Pb it is still Cn the dominant contribution, 
with an intermediate region around 1.5 A~^, where Dn 
is of similar magnitude. Therefore the overall decay of 
Fii (fc, t) is dictated by z^^^ , which being much larger than 
z*^^^ induces a much faster decay; note the smaller scale 
in the time axis for Fii{k,t) in fig.|S| 

It is also interesting to observe that out of the two 
propagating roots, the dominant one for very small k 
is the low frequency one, and therefore in this limit all 
the three partials oscillate with the same frequency, as 
happens within the hydrodynamic model. There is how- 
ever a transition from this behavior to the dominance of 
the high frequency mode for k values around 0.2 
for Li4Pb and around 0.5 for Lio.yMgo.a . This will 
have a direct impact on the presence of low and / or high 
frequency peaks in the partial DSFs. 

In fig. ^1 we show the DSFs for both alloys at a small 
wavevector, a wavevector in the transition region and at a 
larger k. In the case of Li4Pb we observe at the smallest 
k a clear peak in S'ii(fc,Ci;) and 5*12(^,0;) and a shoul- 
der in the 5*22(^,1^), all at the same frequency, akin to 
hydrodynamic sound propagation, even though the pres- 
ence of the diffusive term masks the propagating mode 
and S22ik,Lu) shows no side peak. For larger k we find 
a prominent peak in S'ii(fc,a;) due to the high frequency 
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contribution, while there is no indication of any peak in 
£'22(^,0'). This appearance of high frequency peaks in 
the DSF of the hght component (11) and its absence in 
the DSF of the heavy component (22) has been associated 
with the so cahed "fast sound" phenomenon, which as we 
see appears naturaUy within the viscoelastic model as a 
consequence of the two propagating roots of the model. 
In the case of Lio.yMgo.s the diffusive contribution to 
Sij{k,uj) is very large for small k, and completely covers 
the (low frequency) propagating contribution, so that no 
side peaks appear in the DSFs even for rather small k 
values. Anyhow this contribution is indeed present and 
shows up as a rather weak shoulder at the same frequency 
in all the three Sij{k,uj). For larger k the higher fre- 
quency component becomes more important in 5*11 (fc,w) 
and due to its larger frequency is less covered by the dif- 
fusive contribution, so that for sufficiently large k a peak 
finally develops in 6*11(^,0;), whereas S22ik,uj), whose 
dominant propagating contribution is the low frequency 
one, never develops a side peak. Therefore the situation 
is rather similar to the case of Li4Pb, except that the 
separation between the frequencies of the low and the 
high frequency modes is smaller due to the smaller mass 
of Mg, and that the magnitude of the diffusive mode in 
Lio.rMgo.s definitely obscures the analysis. 

Coming to the BT functions, one can expect that 
Lio.vMgo.s , being a quasi-ideal system, should behave 
similar to the pseudobinary system, i.e., the single- 
particle dynamics should dominate the concentration- 
concentration partials, while the collective dynamics 
should dominate the number- number partials. This is 
indeed what happens, as shown in fig. 1141 where we see 
that FN]y{k,t) in fact decays much faster than the par- 
tials, which is similar to the behavior obtained for the 
pseudobinary system, Fcc{k, t) is practically purely diffu- 
sive, as happens for the SISF of a one-component system, 
and Ff^[c{k,t) is very much smaller than the two other 
partials. Obviously this is also reflected in the dynamic 
structure factors, which are plotted in fig. El SNNik,uj) 
shows clear shoulders or peaks, while Sccik, w) shows nei- 
ther side peaks nor shoulders (note that for the smallest 
k shown the MD FMc{k, t) is too noisy to allow a reliable 
Fourier transform). On the contrary, the case of Li4Pb is 
different, as shown in figs. ITCl and [T7I FNN{k,t) now ex- 
hibits a certain diffusive component and Fcc{k,uj) does 
show oscillations. As a consequence, Scc{k,uj) presents 
well defined side peaks or shoulders, indicative of propa- 
gating concentration modes. It is important to note that 
the magnitude of these peaks is very small, around 10~^, 
compared to the value of Scc{k, uj = 0), which is far out- 
side the range of the graph, being around 10"'^ — 10~^ 
depending on the wavevector. We stress that the repro- 
duction of the different behaviors exhibited by the two 
alloys, even at the level of such small details, represents 
a stringent test which is well fulfilled by the viscoelastic 
model. 

Finally we address the longitudinal current correlation 
functions. In fig. E| we show the AL and the BT func- 



tions for Lio.yMgo.a , and in fig. E| those correspond- 
ing to Li4Pb. Again the case of Lig yMgo a is very sim- 
ilar to pseudobinary Li. For small k, the three AL par- 
tials have a peak at a common frequency. On the other 
hand, Ci2{k,uj) exhibits both a maximum and a mini- 
mum for all fc, whose positions coincide with the max- 
ima (or shoulders at small k) of Cii(fc,Ct') and C22{k,oj) 
respectively. The behavior of the BT currents is also 
similar to pseudobinary Li: CNN(k,uj) only shows one 
peak, and Ccc{k,uj) also shows only one peak, although 
since the decoupling number-concentration is not com- 
plete, a clear feature (a shoulder) does also appear at 
the frequency of the maximum of CNN{k, oj). In the case 
of Li4Pb it is evident that the smallest k shown in the 
graph is already outside the hydrodynamic region. For 
all the k values shown we observe two characteristic fre- 
quencies, a small one for which we have a maximum in 
C22{k,i-u) and a maximum (at low k) or a minimum (at 
larger k) for Ci2{k,iu), and a high one at which there is 
a maximum in Cii{k,uj) and a minimum (at small k) or 
a maximum (at larger k) in Ci2(fc,w). The BT currents, 
on the other hand, always show a common maximum 
at the high characteristic frequency and a maximum in 
Ccc(fc,w), minimum in CiVc(fc,w), at the low one. We 
therefore see that in Li4Pb Ccc{k,uj) has a clear double 
peak structure, each situated close to the low and high 
frequency roots of the model. 



IV. CONCLUSIONS. 

We have presented an extension of Lovesey's viscoelas- 
tic model to binary mixtures that correctly recovers the 
same model for the one-component system in the case of 
pseudobinary mixtures, including both the collective and 
the single-particle properties. 

We have observed in the MD simulations of this kind 
of system a clear appearance of propagating modes in the 
single-particle dynamic properties of the one-component 
system, specifically in the longitudinal self-current corre- 
lation function, which the viscoelastic model reproduces 
correctly. 

In the case of alloys, the appearance of two propagat- 
ing roots in the model leads naturally to the fast sound 
phenomenon, which is explained as the dominance of the 
high frequency propagating contribution in the DSF of 
the light component, appearing near the k values where 
the amplitudes of the two propagating roots coincide. For 
smaller k the amplitude of the high frequency component 
goes to zero and therefore all the DSFs have a compo- 
nent with a peak at a common frequency, namely that of 
the smaller propagating root. Whether this component 
shows up as a peak or not in the DSFs, depends strongly 
on the amplitude (and width) of the diffusive contribu- 
tions, which in some cases (like Lio.yMgo.s ) can cover 
the propagating contribution. 

The viscoelastic model is able to reproduce the dif- 
ferent behaviors of the ISFs of pseudobinary Li and 
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Lio.yMgo.s on one hand and Li4Pb on the other. The ex- 
tremely slow decay of the ISFs for the former two cases 
is dictated by the small value and large amplitude as- 
sociated to the smallest real root of the model, which 
moreover in the case of pseudobinary Li is strictly the 
single-particle contribution. In the case of Li4Pb a com- 
bination of factors leads to a more rapid decay: first, the 
smallest real root is larger than in Lio.yMgo.s ; second, 
the amplitudes associated to this root in the 22 and 12 
partials are much smaller than those in Lio.yMgo.s , due 
to the behavior of the amplitudes which follow the corre- 
sponding partial structure factors; and third, in the case 
of the 11 partial it is the larger real root the dominant 
one for small k, leading therefore to a much faster over- 
all decay. Another important difference is the presence 
in Li4Pb of propagating concentration modes, which are 
absent in Lio.yMgo.s . The presence (or absence) of clear 
side peaks in Scc{k,oj) is neatly reproduced by the vis- 
coelastic model, in spite of its very small magnitudes. 

We conclude by stressing that although quantitative 
differences do appear between the results of the model 
and the simulations, nevertheless the model reproduces 
accurately the main characteristic features of the dy- 
namic properties of the systems studied, including very 
fine details as the peaks of Scc{k, w) and even unexpected 
peaks in the single-particle properties. We believe that 
this warrants a wide application of the model to the 
analysis of data (possibly obtained by any other means) 
on the dynamic properties of alloys, and on the single- 
particle properties of one-component systems, which will 
complement the already usual application of the model 
for the collective dynamics of one-component systems. 
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APPENDIX 

Here, we derive the expressions for the matrices H and 
K(fc) in the case of a two-component system. In Section 
IIIBI these matrices were determined by the integral 



dt exp 



-N(fc,t = 0) 



The exponential is easily performed after diagonaliz- 
ing N(fc,t ^ 0). Denoting by {Ai(fc), A2(fc)} and 
{xi(fc), X2(fc)} its eigenvalues and eigenvectors we can 
write 



where Nd (k) is a diagonal matrix having the eigenvalues 
{Ai(fc), A2(fc)} on the diagonal, whereas P(fc) is a matrix 
having {xi(fc), X2(fc)} as its column vectors. The expo- 
nential is then straightforward and the wanted integral 
is readily evaluated leading to 



M(fc,z = 0) = ^ H P(fc) D(fc) P(fc) ^M(fc,< = 0) 

(A.2) 

where D(fc) is a diagonal matrix with elements 
{^J2/\l{k), ■\/2/A2(fc)} on the diagonal. Now, by using 
eqn. ifT^ . it is obtained that 



F(fc,z = 0) = 

^ [HP(fc)D(fc)P(fc)-iM(fc, i = 0)1 ^ F(fc, t = 0) 

(A.3) 

and by examining the /c — s- oo limit of this equation, the 
elements of the matrix H are obtained. In particular we 
have 



Xi{k oo) 
M{k^oo,t^ 0) 

J ij 

F(A: -> oo,z = 0) 



2k^kBT 



PksT 



nSij (fc, Lo = 0) = 

1/2 



2fc2 ksT 



(A.4) 



leading to H = (2\/2/7r)I. 

Now, we turn to the determination of the matrix K(fc). 
By inserting eqn. ljA.2|) in eqn. (|25|l . and using the pre- 
vious result for H, it is obtained that 



K(fc) = ^ P(fc) D(fc) P(fc) ^ (A.5) 

The explicit expressions for the matrix elements of K(A:) 
are 



2 iVi2(fc)A^2i(fc)Ay'(fc) +/3^(fc)A2^^(fc) 

A iVi2(fc)iV2l(fc)+/3'(fc) 

2 NMf3{k)[Xy\k) - Xl/\k)] 
A NMk)N2i{k) + I3^k) 

2 N2i{k)m[X'/\k)-Xl^^(k)] 
A Ni2{k)N2i{k) + (3^k) 
2 Ni2{k)N2i{k)Xl^\k) + p\k)X\^\k) 



1/2, 



Kii{k) = 

Ki2{k) = 

K2iik) = 

K22{k) = 



iVl2(fc)iV2l(fc)+/32(fc) 



N(fc, i = 0) = P(fc) ND(fc) P(fc)" 



(A.l) 



(A.6) 
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where the Nij{k) are the matrix elements of N(fc,<: = 0), 
the eigenvalues are given by 



Hk) = ^[Nn{k) + N22{k)] 



Nii{k) - N22{k) 



Ai(fc) = 7^[Nii{k) + N22{k)] 



+ Ni2ik)N2i{k) 



+ 



Nn{k)-N22{k) 



+ m2{k)N2i{k) 



1/2 



1/2 



(A.7) 



and p{k) = Xi{k) - Nu{k). 
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FIG. 1: Static partial structure factors for pscudobinary Li. 
Full line and circles: Sii{k) = 522 (fc), dashed lino and trian- 
gles: Si2{k), dotted line: one-component system S(k). The 
lines are theoretical results and the symbols denote MD sim- 
ulation results. 
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FIG. 3: Roots of pseudobinary Li in the viscoelastic model. 
Full lines: Jf' and F^^^ dashed Unes: ujf and F^^^ dotted 
line: dash-dotted line: z'^'^^ . 
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FIG. 2: Partial ISFs for pseudobinary Li. Full lines (up- 
per panel): Fii{k,t) = F22{k,t), dashed line (mid panel): 
F\2{k,t), dotted line (lower panel): FNN(k,t), which coin- 
cides with the one-component system F{k,t). The symbols 
are the corresponding MD results. Left column: k = 0.23 
A"'^, mid column: k = 1.23 A~^, right column: k = 2.51 
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FIG. 4: Amplitudes of the different modes of pseudobinary Li 

in the viscoelastic model. Full lines: A\\, dashed lines: Bii, 
dotted line: Cn, dash-dotted line: D\\. 
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FIG. 5: Partial currents for pscudobiriary Li. Lines are the- 
oretical results and symbols denote simulation data. Full 
line and circles (upper panel): Cii{k,Lo) = C22{k,uj), dashed 
line and triangles (upper panel): Ci2(fe, w), dotted line and 
circles (lower panel): CArjv(fe,w), which coincides with the 
onc-componcnt system C{k,u)), dash-dotted line and trian- 
gles (lower panel): scaled Ccc{k,uj), which coincides with the 
onc-componcnt system Cs{k,uj). Left column: k — 0.23 A""', 
mid column: k = 1.23 A~^, right column: k = 2.18 A~^. 




FIG. 6: Static partial structure factors for Lio.rMgo.a and 
Li4Pb. Full lines and circles: S'n (fc), dotted lines and squares: 
522 (fc), dashed lines and triangles: Si2(fc). The lines are theo- 
retical results and the symbols denote MD simulation results. 
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FIG. 7: Partial ISFs for Lio.yMgo.a. Linos arc theoretical 
r(>sults and symbols are sinnilation data. Full lines (upper 
panel): Fii{k,t), dEished line (mid panel): F\2(k,t), dotted 
line (lower panel): F22{k,t). Left column: k = 0.28 
(theory), 0.26 (simulation); mid column: k = 1.25 
(theory), 1.23 A"'^ (simulation); right column: k = 1.88 A^^ 
(theory and simulation). 




FIG. 8: Partial ISFs for Li4Pb. The meaning of the lines and 
symbols are the same as in the previous figure. Left column: 
k = 0.31 A~^ (theory), 0.27 A~^ (simulation); mid column: 
k = 1.30 A^^ (theory and simulation); right column: k — 1.76 
A~^ (theory), 1.74 A~^ (simulation). Note the different scales 
in the time axis for Fn {k, t). 
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FIG. 9: Roots of Lio.rMgo.s (thick lines) and Li4Pb (thin 
Unes) in the viscoclastic model. Full lines: ajs^' and r'^\ 
dashed lines: and F*^', dotted line: a'^-*, dash-dotted 

line: z^''\ 



!// 

'^^ 


0.04 


^ y / 
/■-' , */ - 




/ ,->■"' 


0.02 / 




// / \ 


0.5 1 1 


• — J-"-^* 


L\>', .L^ :-.-t-£^ ^ :.i»„ 



12 3 4 



FIG. 10: Amplitudes of the different modes of F\\{k,€) for 
Lio.rMgo.a (thick linos) and Li4Pb (thin lines) in the viscoelas- 
tic model. Full hncs: An, dashed lines: Bn, dotted line: Cn, 
dash-dotted line: D\\. 
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FIG. 11: Amplitudes of the different modes of Fi2{k,t) for 
Lio.rMgo.a (thick lines) and Li4Pb (thin linos) in the viscoelas- 
tic model. Full lines: A 12, dashed lines: B12, dotted line: C12, 
dash-dotted line: Di2- 
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FIG. 12: Amplitudes of the different modes of F22(fe,t) for 
Lio.rMgo.a (thick hnes) and Li4Pb (thin Hnes) in the viscoelas- 
tic model. FulUincs: A22, dashed lines: B22, dotted line: C22, 
dash-dotted line: 7322 • 




FIG. 13: Dynamic structure factors for Lio.rMgo.a (upper 
pannel) and Li4Pb (lower panncl) at wavevectors shown in 
the graphs. Full lines: (fc, w), dashed lines: Si2(fe, w), dot- 
ted line: ^22 (A;, w). The thin dash-dotted line denotes the low 
frequency contribution to Sw (k, uj) and the thin full line is the 
corresponding high frequency contribution, which is negligible 
for the lowest k. Symbols denote simulation data at wavevec- 
tors 0.53, and 0.65 A^^ for Li4Pb. For those graphs S22{k,oj) 
have been shifted upwards by 0.0001 to improve visibility. 
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FIG. 14: Bhatia- Thornton ISFs for Lio.rMgo.s . The k values 
for the different columns are the same as in fig. Q Note in 
some cases the different scales in the time axis for Fjvjv(fe,t). 
Lines are theoretical results and symbols are simulation data. 




FIG. 15: Bhatia-Thornton DSFs for Lio.yMgo.s . The k values 
for the different columns are the same as in fig. |7| Lines are 
theoretical results and symbols are simulation data. 




FIG. 16: Bhatia- Thornton ISFs for Li4Pb. The k values for 
the different columns are the same as in fig. |H] Lines are 
theoretical results and symbols are simulation data. Note in 
some cases the different scales in the time axis for _Fjvjv(fc, t). 




FIG. 17: Bhatia- Thornton DSFs for Li4Pb. The k values 
for the different columns are the same as in fig. |H| Lines are 
theoretical results and symbols are simulation data. 
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FIG. 18; Longitudinal current correlation functions for 
Lio.7Mgo.3 . Lines are theoretical results and symbols are sim- 
ulation data. Upper pannel: Ashcroft-Langreth partials; full 
line and circles: Cii(fc,aj), dashed line and squares: Ci2{k,ui), 
dotted line and triangles: C22{k,uj) Lower pannel: Bhatia- 
Thornton partials; full line and circles: CNN{k,uj), dashed 
line and squares: scaled Ccc{k,uj), dotted line and triangles: 
scaled Civc(fc, w). The k values for the different columns are 
the same as in fig. Q 
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FIG. 19: Longitudinal current correlation functions for Li4Pb. 
Lines are theoretical results and symbols are simulation data. 
Upper pannel: Ashcroft-Langreth partials. Lower pannel: 
Bhatia-Thornton partials. The meaning of the lines and sym- 
bols is the same as in the previous figure. The k values for 
the different columns are the same as in fig. |2l 



